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ABSTRACT 

A simple radiative transfer model of the particle layer found at 
85 km over the summer poles is presented. The effects of the layer on 
the global radiative temperature, the polar region temperature, and the 
greenhouse effect are discussed. The estimated magnitude of the global 
radiative temperature change is 3. 5 x 10 K to 2. 2 x 10 K, depending on 
the value of the imaginary part of the particle index of refraction. The 
layer is shown to have a possible secondary influence on the temperature 
of the polar region while the contribution which the layer makes to the 
greenhouse effect is shown to be negligible. The imaginary part of the 
particle index of refraction is shown to be important in determining the 
attenuation properties of the layer. 



CHAPTER I 


INTRODUCTION 

Attention has recently been focused upon a layer of particulate 
matter found at an altitude of about 85 km. As reported by Donahue et al. , 
(197E), the layer, which is believed to consist of ice crystals, is found 
over the geographic poles during the local summer months- Noctilucent 
clouds, also found in the same altitude region, are thought to be extensions 
of the layer - 

Although little is presently known about the layer, it could play a 
significant role in processes of the region- As shown in a recent paper by 
Olivero (1974) surface reactions on these upper atmospheric aerosols could 
affect the local photochemical balance. The radiative transfer properties 
of the layer could help to explain the mystery of the extreme temperatures 
of the high latitude mesopause and influence polar climatology in a 
non-negligible manner . 

This report examines the radiative transfer properties of the layer 
with regards to global climatic effects- The model developed represents 
one part of a larger study which will encompass both the chemistry and 
dynamics of the region. As expected, the particle layer was found to have 
a negligible effect on the global effective radiative temperature of the 

Earth, the estimated magnitude of the change ranging from about 3. 5 x 

-3 -2 

10 K to 2.2 X 10 K. The layer is shown to have a possible secondary 
influence on the temperature of the polar region. Also, the contribution 
which the layer makes to the greenhouse effect is shown to be negligible. 
The extent to which the attenuation properties of the layer are dependent 
upon the imaginary part of the particle index of refraction is also shown. 
Finally, the composition and the shapes (geometry) of the particles 
making up the layer are discussed. 



CHAPTER n 


ATTENUATION CHARACTERISTICS 
2. 1 Basic Considerations 

Radiation interacting with a particle can either be scattered by the 
particle, at an angle© from the direction of propagation,' or it can be 
absorbed by the particle. Here we assume the particles to be spherical 
and to be evenly distributed within a layer of thickness 6z. The layer is 
sufficiently tenuous to allow multiple scattering to be neglected. 

The radiation incident upon the layer, I^, can be reduced by 
scattering and absorption within the layer. The radiation directly trans- 
mitted by the layer, is the incident radiation less that which has been 
scattered away and absorbed by the particles. 

The scattered radiation can be broken down into forward and back- 
ward components. The backward component, I|^g» represents that radiation 
reflected back to space. The forward component, 1^^, represents that 
passed through the layer but in directions other than the direction of 

propagation of the incident radiation. At a large enough distance below 

» 

the layer, Az » 6z, the total intensity one would measure, I^ , is equal 
to the sum of that transmitted and that scattered in the forward direction. 


I 



h + h 

t fs 


{ 1 ) 


Calling the intensity of the radiation lost to the particles via 

absorption we can express the radiation budget as 


I 

o 


+ L + I - 
t bs abs 


( 2 ) 
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2. 2 Determination of Backscattering, Absorption, and Extinction 
Cross Sections 

Let the total energy scattered in all directions by a single particle 
be equal to the energy of the incident wave falling on the area cr^. The 
total energy absorbed is equal to the energy incident on the area 
The total energy removed. from a pencil beam due to the presence of a 
particle is then equal to that incident on the area cr^ where 


The areas cr , cr , and cr , are, respectively, the cross sections for 
e s 2 tbs 

extinction, scattering, and absprption. Recalling that scattering occurs in 
both forward and backward directions we can write the scattering cross 
section as 


tr = (T, + cr,; . 

s bs fs 

Since, for the layer configuration used in. this model, the radiation lost 
is due to backscattering and absorption the effective extinction cross 

T 

section, cr^ , is given by 


1 

& 

e 


'’^bs ^ '^abs * 


(3) 


For a scatterer with a complex index of refraction m, given by 
m = n^ - i n^ , 


the scattering cross section may, in general, be expressed as (Van de 
Hulst, 1957) 




H 

0 (p 


_L 


F (X, 0 , <p ) sin 9 d<^d9 


(4) 
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where K is the wavenumber 


K = 


2tt 

\ ’ 


( 5 ) 


X. being the wavelength of the incident radiation, and F (X, 0, is a. 
complicated function, defining the intensity of the scattered radiation from 
a particle, described by a size parameter X for an arbitrary scattering 
angle 0 and azimuth angle <p . The size parameter is given by 


X = 



( 6 ) 


r being the radius of the particle. The general expression for or given 

s 

by (4) can, for various sizes of X, be replaced by simpler expressions. 
For example, when X « 1, or when the particle is much smaller than 
the wavelength of the incident radiation the Rayleigh scattering theory is 

applicable. In this case the scattering cross section is found to be 

-4 

proportional to X . (See Equation (17).) 

When the particle is large compared to the wavelength, X » 1, 
the general theory according to Mie is employed. In this case F (X, 0, 
<f?) may be expressed as 


2 2 

F (X, 0, <p) = i^tX, 0) cos <p + ij (X, 0) sin <p 


(7) 


where ij^ (X, 0) and ^2 ©) ^^6 the intensities of the scattered radiation 

per solid angle perpendicular and parallel, respectively, to the plane of 
scattering. Combining Equations (4) and (7) we get 


(M 






iy (X, 0) cos (p + i, (X, 0) sin cp] 


0 <p 


sin 0 d(p d0 . 


( 8 ) 
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Assuming the radiation to be linearly polarized and substituting Equation 
(5) for K, we can express ( 8 ) as 

(X) = ^ y [i^ (X, e) + i 2 (x, 0)] sin e de . (9) 

e 

The intensity functions ij^ (X, 0) and i 2 (X, 0 ) are complicated expressions 

involving the Riccati-Bessel functions and their derivatives. When the 

complete expressions for ij^ (X, 0) and i 2 (X, 0 ) are substituted into (9) 

-2 

the Mie scattering cross section is found to be proportional to X 

The absorption cross section is more difficult to represent. In 
Van de Hulst's (1957) treatment an absorption cross section is developed 
but the phenomenon represented by the cross section is backscattering 
plus pure absorption. This definition has been used in aerosol models 
(Ensor et al. , 1971) and is a feature of a popular Mie scattering program 
(Dave, 1968 ) but this definition tends to obscure what is meant by absorption 
by a particle. The absorption properties of a particle are best treated 
theoretically from a quantum mechanical standpoint or experimentally 
(e.g., Lin, Baker, and Charlson, 1973). The absorption cross section 
can be approximated for small X as (Lin, Baker, and Charlson, 1973) 

' '”1- "z' 


where n^ and n^ are, respectively, the real and imaginary components of 
the index of refraction. For values of X < 10, Equation (10) can be 
approximated as 


abs ' ' 


^23 
27 t r 


n^ f (nj^) 


( 11 ) 


For a real component of index of refraction n^ = 1.33, f (nj^) is assigned 
the value of 1. Therefore, the absorption cross section is given by 



( 12 ) 




-23 
2tt r 


n. 


In the present study the imaginary component of the index of refraction 
was varied from 0 to 0. 9. 

2. 3 Optical Depth Determinations 

Consider the layer geometry shown in Figure 1. As shown by 
Equation (2) the incident radiation is equal to the sum of that transmitted, 
absorbed, and reflected by the layer. Employing Beer's law the radiation 
transmitted by the layer of thickness 6z is given by 

z) = Iq(M exp [-T (X, z)] (13) 

where t(X, z) is the optical depth. The optical depth for vertically 
incident radiation is given by 

z + 6z 

t(X, z) = y ^ %i (14) 

z i 

I 

where <r^. (X) is the extinction cross section of constituents i, given by 
Equation (3), m (z) is the number density of constituent i, and the 
summation is performed over all of the constituents making up the layer. 
It is assumed here that the layer consists of both particulate matter and 
the molecular constituents making up the atmosphere. For the high 
latitude particle layer under consideration the incident radiation is not 
normal to the layer, therefore, the zenith aingle, x> of th® radiation must 
be considered when determining the amount of transmitted radiation. As 
X increases, the path length of the radiation increases and the amount of 
transmitted radiation decreases. Therefore, the optical depth becomes 



I 



Figure 1. Layer geometry 
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z + 6z 

t{\ , z, X) = sec X I I dz . (15) 

z i 

As stated previously, the layer has been found to extend over the summer 
geographic poles. For this work, it is assumed that the layer extends 
from 75° latitude to the poles. Furthermore, it is assumed that the 
number density of the particulate matter is uniform throughout the layer. 
With these assumptions and approximating sec x ® sec x> (15) can be 
written as 

z + 6z 

t(X, z) = sec X y ^ {^) n^(z) dz (16) 

z i 


where x Is the mean zenith angle for the layer. Table 1 lists the 
properties adopted for the layer. 

The adopted atmospheric model consists of gaseous constituents, 
notably O, O 2 , and N 2 , and the particles making up the layer. It is 
assumed that for wavelengths above about 0.2|j, the gaseous (molecular) 
constituents interact with solar radiation only by scattering. (Note that 
at these high zenith angles little radiation at wavelengths below 0, Z[x 
reaches the mesosphere.) The scattering cross section for these 
constituents, describable by the' Rayleigh scattering theory is (Van de 
Hulst, 1957) 





X n 


m 



6 + 36 
6-76 


(17) 


where m is the index of refraction of the atmosphere, assumed real, 

a, 

is the total number density of the molecular constituents, and 6 is the 
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Table 1. Adopted Properties of the Particle Layer 


Altitude, z 

85 km 

Thickness, 6z 

3 km 

Spatial 

Uniform from 75° latitude 
to the poles 

Number density of molecular 
constituents, n 

m 

14 -3 

1. 532 X 10 cm 
(CIRA, 1965) 

Aerosol number density, n^ 

33 cm”^ 

Aerosol particle radius 
(monodispersed), r 

0. 13p 

Mean zenith angle, X 

79.7° 


Range of particle index of 
refraction, m 


1. 33 - 01 to 1, 33 - 0. 9i 
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depolarization factor, assumed negligible. The index of refraction is, 
in general, given by (Bean, 1962) 


, . 7.76x 10"^P . 7.72xl0‘^e . 0.375 e 

^ T T + "^2 


(18) 


where P is the pressure in millibars, T is the temperature in degrees 
Kelvin, and e is the partial pressure of water vapor in millibars. 
Assuming negligible amounts of water vapor we can approximate (18) as 


m “ 1 + 
a 


7.76 X 10" P 


(19) 


>-9 


For this work m » 1, m - 1 = 1.66 x 10 . In Figure 2 the Rayleigh 

a a 

cross section is plotted against wavelength. The optical depth for the 
molecular constituents may be expressed as 


z) = sec X 


O' (X) 
om ray 


( 20 ) 


where H is the mean scale height of the atmosphere and n is the 

a or Qj^ 

number density of the molecular constituents at the base of the layer. 
The optical depth expressed by (20) is that for a column of air extending 
from the base of the aerosol layer to infinity. 

For the aerosol components the optical depth is given by 


(X, z) = sec X tr' (X) n^ 6z (21) 

JLi C X-r 

where n_ is the average number density of the particles in the layer* 

XJ 

From Equation (3) we can write (21) as 


(X, z) = sec X (X) + (X)] n^^^ 6z (22) 

where (X) is the integrated backward scattering cross section referenced 
to the layer orientation as shown in Figure 1. For this orientation (X) 
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~ ^bs Figures 3 and 4 are, respectively, the backscattering and 

absorption cross sections plotted against wavelength for various values 

of n 2 » It is evident that as the imaginary part of the index of refraction, 

U2t increases the absorption cross section begins to dominate. This 

result has been confirmed by other workers (Kattawar and Plass, 1967; 

Plass, 1966; and Ensor et al. , 1971). As one might expect, the particle 

layer optical depth was found to be greater than that of the molecular 

constituents above the 85 kh-j level, the ratio /t being on the order of 

Ju/ rn 

10^ to 10^. 

In Figure 5 the layer vertical optical depth is plotted against 

wavelength for various values of n 2 « In Figure 6, the layer vertical 

optical depth is plotted against n 2 for a wavelength of 0. 55|x. It is again 

worth noting that the backscatter cross sections used in determining (22) 

are integrated over all backscatter angles. The results reported by 

Donahue et al. , (1972) were obtained from the OGO-6 photometer which 

measured airglow emissions of sodium (5893 ^) and the atomic oxygen 

green line (5577 ^), (Donahue et al. , 1973). Their vertical optical depth 
-5-5 

values of 4 X 10 to 7 x 10 were determined for a scattering angle of 
80°, particle radius of 0. 13(j., and real index of refraction of 1.33 
(Guenther, personal communication). Using these conditions and a 

_4 

wavelength of 0. 5577|j., the vertical optical depth from (22) is 1. 5 x 10 
From (22) and knowing that Tj^ » we can write (13) as 

Ij (z, X) ~ I^(X) exp {- sec X [ (X) + (X)] n^ 6z}. 

(23) 

By dividing both sides of (23) by Iq(X-) one 'gets an expression for the 
transmittance of the layer, T^, the fraction of the incident radiation 
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Figure 3. Integrated backscatter 
wavelength 


= 1.33 -0.9 i 


I 

10.0 

cross section vs. 


cTq (cm2) 
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Figure 4* Absorption cross section vs* wavelength 
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that is transmitted by the layer, 

yz, X) 

Tr(z, X) = -wxF" {-sec x 

° (24) 

In Figure 7,1- (z, X), the effective layer extinction, is plotted 

against for X = 0. 55|j.. The first feature that is evident is that the 
layer is essentially transparent to the incoming radiation. For longer 
wavelengths, X> l.Op, the transparency is more striking. The other 
feature is that as increases T^ decreases, due to the increasing 
importance of absorption within the layer. 

One may also determine the layer albedo, , the fraction of 
incident radiation reflected back to space, 

{z, X) = (z, X) / I^(X), 

where (z, X) is the reflected radiation, given by 

1^^ (z, X) = (1 - exp {- sec X (X) n^^ 6z}). (25) 

—X 

Employing the approximation that e = 1 *■ x we can express the 
albedo as 

A^{z, X) « secx (T^g (X) n^^ 6z . (26) 

In Figure 8 the layer albedo is plotted against for \ = 0. 55|jl and x “ 

79.7^» Again, the importance of n^^ and hence absorption within the 

layer is evident. The vertical optical depths of Donahue et al. , (197Z) 

-4 -4 

correspond to albedo values of 2. Z x 10 to 3. 9 x 10 when the zenith 

angle dependence is included. For noctilucent clouds, Sharanov (I960) 

-5 -5 

has calculated apparent albedo values of 2.3 x 10 to 4.7 x 10 
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Figure 7. Effective layer extinction vs. 


- 19 



Figure 8. Layer albedo vs. 
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Employing a scattering angle of 80*^ and a wavelength of 0. 5577jj. the 

_4 

albedo calculated via (26) is 8.4 x 10 



CHAPTER III 


ENERGY CONSIDERATIONS 


3. 1 Global 

Consider now the balance between the radiant energy absorbed by 
the Earth and that emitted by the Earth. This may be expressed as 

(1 - A) S^ = 4tr T^ (27) 

where A is the planetary albedo, S^ the incident solar flux striking the 
atmosphere, tr the Stefan-Boltzmann constant, and T the mean radiative 
temperature. It should be realized that the mean radiative temperature 
is not the surface tempcirature of the Earth but that which would be 
obtained from radiometric observations from outside the atmosphere. 
The solar flux values used in this work were taken from the Handbook 
of Geophysics and Space Environment (1965) and are listed in Table 2. 

A planetary albedo of 0. 33 was taken from Pollack and Toon (1974). 

Let us assume that Equation (27) represents the case where the 
aerosol layer is absent. Assuming, now, the presence of the layer we 
can write for the energy balance 

(1 - A) S = 4 ir (T')^ (28) 

where S, the mean flux perturbed by the presence of the aerosol layer, 
is given by 

S = aS. + bS (29) 

J— f o 

where a is the fractional area of the Earth covered by the particle layer, 
Sj^ is the solar flux transmitted by the layer given by (22), b is the 
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Table 2. Solar Fluxes Incident on the Top of the Atmosphere 
(Handbook of Geophysics and Space Environments, 1965) 


\ (|j,) 


Flux (watts m 



0. 27 
0. 28 
0. 30 
0.32 
0.34 
0. 36 
0.38 
0.40 
0.44 
0. 50 
0. 55 
0. 65 
0. 65 

0. 70 
0. 80 
0. 90 

1 . 02 

1. 25 
1. 60 

2. 50 

3. 50 

4. 25 


2. 30 
2. 60 
14. 91 
17. 97 
22 . 62 

24. 96 

25. 14 
30. 96 

129. 28 
100. 50 
97.06 
91,53 
81.31 
105. 18 
113. 09 
89.79 
104. 59 
136. 77 
113.97 
53.73 
15.34 
3.45 


* Integrated over adjacent wavelength intervals, A\, 
indicated by interval centers 
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fractional area of the Earth not covered by the layer, and T' is the mean 
radiative temperature in the presence of the layer. Knowing that b = 1 - a, 
we can write the temperature difference, AT, as 


AT 




(30) 


The percent change ifi solar fl\uc due to the layer. 



“3 

is shown in Figure 9 plotted against n. 2 ^, ranging from 5. 6 x 10 %to 

- 2 -3 

3, 5 X 10 %, The global temperature difference ranged from 3. 5 x lO'^^K 

_2 

to 2.2 X 10 K, again depending on shown in Figure 10. The 
temperature difference averaged globally represents a negligible cooling 
of the Earth. Again, it is mentioned that the temperature changes reflect 
changes in the mean radiative temperature of the Earth and not the surface 
temperature. In order to make surface calculations the energy of the 
atmospheric circulation, condensation, and evaporation must be considered. 
As mentioned by Rasool and Schneider (1971) surface temperature changes 
on the order of 3. 5 K, if sustained over a period of several years, could 
be enough to trigger an ice age. The significance of such particle layers 
on global climatology is beyond the scope of the present investigation. 


B. Polar 

The radiative balance of the polar region can also be examined 
using a form of (27). In this case the radiative energy balance is given by 


a 

c 


P o 


4 

a 4 0 - T 
s 


(31) 


where a^ is the fraction of the projected area poleward of 75° latitude. 



AS/So(%) 




AT V8 
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A is the albedo of the polar region, taken to be 0.75 {Atwater, 1970), 

P ■ 

and a is the fraction of the surface area poleward of 75° latitude. Again, 
s 

a temperature difference between the layer and non-layer cases can be 
found, given by 


AT 


r a (1 - A ) 1 
c ' P . 

I 

4 

" 1 

1 


C 4 

S 4 

L a 4 (T 
L s 


o 

O 

L J 


(32) 


In Figure 11, AT is plotted against for a zenith angle of 56,4^, 

-2 

corresponding to the summer solstice, and ranges from 6. 09 x 10 K to 
0. 533 K. The low n^ values represent insignificant temperature 
differences but the upper extreme values are possibly non-negligible. 

The layer is a present feature of the earth-atmosphere system and 
represents a cooling term in the energy budgets If the layer was removed 
the polar region would, in this steady state model, experience a warning 
up to a few tenths of a degree, depending on n^p due to the added radiation 
reaching the surface* This warming alone most likely would not be enough 
to alter the climatology of the region, but as Oort (1974) has pointed out 
it may have a secondary influence modifying the temperature anamoUes 
brought about by variations of the energy influx from middle latitudes or 
the Earth^s surface* 

C* Contributions to the Greenhouse Effect 

The radiation emitted by the Earth undergoes attenuation by the 
layer also* Radiation reflected back to the Earth from the layer provides 
another heat source for the overall energy budget of the Earth, contributing 
to the so-called greenhouse effect* As shown by Figure 3, the backscatter 
cross section is extremely small for longer wavelengths, approaching the 
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Figure 11. At vs. (polar average) 
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-15 2 

Rayleigh regime. For a wavelength of 1 1, Oji, 

m = 1. 33 - 01, which corresponds to a vertical optical depth of - 
“8 

10~ . Approximating the transmittance as 

Tj,(z, \) - 1 - 

it can be seen that the layer is essentially transparent to the Earth's 
outgoing radiation. The heating, resxilting from the reflected radiation 
from the scattering layer is on the order of 10 "^K averaged globally. 

_5 

For the polar region the heating is of the order of 10 K. 



CHAPTER IV 


COMPOSITION AND SHAPE OF THE PARTICLES 

The compaaition of the particles in the layer is crucial to under- 
standing the attenuation properties of the aerosol layer. Unfortunately, 
little is known about the precise make-up of the particles. Polarization 
measurements of noctilucent clouds, which are believed to be equator- 
ward extensions of the layer, .indicate that noctilucent cloud particles 
are non -metallic with a real refractive index between 1.2 - 1.6 (Witt, 
i 960 a, b; Villmann, 1962). As mentioned earlier, the imaginary part of 
the index of refraction is more difficult to establish. 

Rocket samples of noctilucent cloud particles revealed the 
presence of nickel and iron in the particles, suggesting a meteoric 
origin (Witt et al. , 1964). Also, in one series of rocket samples the 
particles were found to have been coated with a volatile substance, at 
the time of collection, believed to be ice, (Hemenway et al. , 1964; 
Skrivanek and Soberman, 1964). The findings of Donahue et al. , (1972) 
involving the particle layer were based on the assumption of spherical 
ice particles of radius 0, 13|ji ahd real refractive index of 1, 33. 

The value of 1. 33 for the real part of the particle index of 
refraction was chosen for two reasons. First, use of n^ = 1.33 allowed 
the results from the model to be compared with the results of Donahue 
et al, , (197 2). Second, tabular results for the scattering functions of 
particles with = 1.33 are known, thereby permitting comparison of 
the results in the model. 

The shapes of the particles are also poorly understood. Classically, 
the particles have been assumed spherical in order to simplify the scatter- 
ing problem. Rocket samples reveal, though, that the particles need not 
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be spherical (Farlow et al. , 1970 and Frank et al. , 1970). Their results 
show that the particles may be rodlike, irregular, or have fluffy, dendritic 
shapes in clusters. 

The shapes of particles would greatly affect the scattering and 
polarization results and the area-to-mass ratio data for meteorites. As 
mentioned before, theoretical results of scattering and polarization 
assume spherical particles for simplicity, the reason being that deter- 
mining scattering from non-spherical particles is an extremely difficult 
task. Also, it is extremely difficult to determine particle shape from 
experimental data (e. g. , lidar), -which then makes it difficult for one to 
make estimates of area-to-mass ratio. 



CHAPTER V 


CONCLUSIONS AND RECOMMENDATIONS 

A simple radiative transfer model has been developed to examine 
the radiation attenuation properties of the particle layer found at 85 km 
over the summer geographic poles. As expected, the layer was found to 
have a negligible effect on the global radiative temperature of the Earth. 

The layer, though, does have a possible secondary Influence on the temper- 
ature of the polar region. The contribution to the greenhouse effect was 
also found to be minimal. In all of these findings the imaginary part of 
the particle index of refraction was shown to be important in determining 
the attenuation properties of the layer. 

The model assumed a monodispersed distribution for the particles 
within the aerosol layer. From polarization measurements (Witt, 1960 a, 
b; and Villmann, 1962) and rocket samples (e, g. , Hemenway et al. , 1964 b) 
of noctilucent clouds it is known that the particles have a size distribution 
that can be approximated by a Junge distribution 

N{d) = Ad"^, 

where N(d) is the number density, d is the particle diameter, A is a 
constant, and 3 < p < 4^ The attenuation properties of the layer assuming 
a polydispersed distribution need to be examined^ 

The greatest source of uncertainty in the present model lies in char- 
acterizing the absorption cross section^ Estimates of the imaginary part of 
the index of refraction, must be considered uncertain to an order of 

magnitudeo In addition, the expression adopted for the absorption cross 

section itself, while probably quite adequate for small values of (of order 
-2 

10 ), is likely to be uncertain to factors of the order of 2 for the larger 

values of used in the present work„ Much still needs to be done to firmly 
establish the absorption properties of upper atmospheric aerosols. 
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